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ABSTRACT 



Aims. The paper presents the timing and spectral analysis of several observations of the Crab pulsar performed with INTEGRAL 
in the energy range 3-500 keV. 

Methods. All these observations, when summed together provide a high statistics data set which can be used for accurate phase 
resolved spectroscopy. A detailed study of the pulsed emission at different phase intervals is performed. 

Results. The spectral distribution changes with phase showing a characteristic reverse S shape of the photon index. Moreover 
the spectrum softens with energy, in each phase interval, and this behavior is adequately modeled over the whole energy range 
3-500 keV with a single curved law with a slope variable with Log(E), confirming the BeppoSAX results on the curvature of 
the pulsed emission. The bending parameter of the log-parabolic model is compatible with a single value of 0.14±0.02 over all 
phase intervals. 

Conclusions. Results are discussed within the three-dimensional outer gap model. 
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1. Introduction 

The Crab pulsar (PSR B0531+21) can be observed in al- 
most every energy band of the electromagnetic spectrum. 
Its pulse profile is characterized by a double peak struc- 
ture with a phase separation of 0.4 that is almost aligned 
in absolute phase over all wavelengths ijRots et al J 12004 
iKuiper et aLll2003t iTennant et aT1l200l[) . 

In the X-ray range, the relative intensity, height and 
width of the two peaks vary with energy: the first peak 
(PI), dominant at low X-ray energies, becomes smaller 
than the second one (P2). Moreover, an enhancement with 
energy of the bridge between these p eaks, usually called 
Interpeak (Ip) , is also well evident l|Mineo et alJ Il997l 
see also Fig.l). At energy above 1 MeV, the morphology 
changes abruptly: the first pulse becomes again dominant 
over the second one and the bridge emission loses signifi- 
cance; the pulse profile above 30 MeV is similar to th e one 
observed at optical wavelengths l|Kuiper et al.ll200l|h 
A first detailed study of th e phase-resolv e d X-r ay spec- 
tra has been performed by IPravdo et all l)l997lh in the 
5-200 keV energy interval, based on RXTE (PCA and 
HEXTE) data. Their main result was a variation of the 
photon index as function of the pulse phase with a reverse 
S shape: the spectrum softens starting from the leading 
edge of the first peak where it reaches the maximum value, 
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it hardens in the interpeak and softens again in the sec- 
ond peak. The S sha pe spectral v ariati on with phase has 
been reporte d bv iMassaro et all ll2000h from BeppoSAX 
data and bv IWeisskopf et alJ l(2004l) from Chandra data, 
even if with lower statistical significance, confirming the 
symmetric evolution of the spectral index around the first 
peak. However, the softening of the PI core r espect to the 
leading edge has been recently questioned by Vivckanand 
(2002) in a new analysis of RXTE data but performed 
over the smaller energy range of 5-60 keV. 
Significant X-ray emission from the pulsar in the off-pulse 
interv al (phase 0.5-0.9) was discovered bv lTennant et alJ 
il200ll) with Chandra observations, however, the spectral 
index measured in t his phase interval suff ers of large sta- 
tistical uncertainty l)Weisskopf et alJl2004J) . 

BeppoSAX observations of the Crab pulsar showed 
that the photon indices of the pulsed emission sig- 
nificantly increase with energy maintaining the same 
S shape behavior over the 0.1-300 keV energy range 
l|Ma,ssaro fc CusumanrjEomt iZhang fc ChendkOO^. The 
spectral index variation has been modeled using a sin- 
gle curved power law with a slope variable with Log(E) 
l|Massaro et al.lEoO(ij) . Moreover, applying this model to 
three wide phase intervals, the first peak, the Interpeak 
and the second peak, a single value of ~ 0.15 for the cur- 
vature parameter has been measured in the three intervals 
( Massaro et~al1l200ll) . 
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50-100 keV (ISGRI) 



100-200 keV (SPI) 




200-360 keV (PICSIT) 





Fig. 1. Crab phase histograms in six energy bands in absolute phase (the main radio pulse at phase 0.0). The light 
curves have different phase resolution according to the statistics avai lable and to the tim e resolution of the instrument. 
PI, P2 and Ip phase intervals according to the definition given in iMineo et al ] lll997l) are indicated in the top left 
panel. 



er et al.l <j200l|) presented a coherent high-energy pic- 
ture of the Crab pulsar from 0.1 keV up to 10 GeV by 
using the high energy 7-ray data from the CGRO satellite 
together with data obtained at X-ray energies from sev- 
eral observatories. The authors model the 0.1 keV-10 GeV 
pulsed emission in 7 narrow phase slices with a composite 
model: a power law present in the phase intervals of the 
two main pulses, a curved spectral component required in 
the same phase intervals and second broader curved spec- 
tral component representing mainly the bridge emission. 

X-ray observations of Crab pulsar performed with a 
balloon born experiment report the detection of an emis- 



sion line a t 440 keV with a flux of (0.86±0.33)xl0 

,-2 „-l 



P h 

cm~" s~" l(Massaro et al.lll99lh . lUlmer et alJ 1^)94), us- 
ing OSSE data, did not detected this line but derived a 3 
a upper limit compatible with its presence. 



Results on the INTEGRAL observation s of the Crab 
pulsar have already been presented by iKuiper et alJ 
ll2003h . that studied the instrument absolute timing ac- 
curacy and bv lBrandt et al.1 l|2003h who reports results on 
the 3-37 keV energy range with the X-ray monitor JEM-X. 

In this paper, we present the timing and spectral 
analysis of several observations of the Crab pulsar per- 
formed with SPI, JEM-X, IBIS/ISGRI and IBIS/PICsIT 



T. Mineo et al.: INTEGRAL observation of the Crab pulsar 
Table 1. Observation log for the data used in this analysis 
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Revolution 


Start-Stop time 




Exposure (ks) 






(MJD) 


JEM-X 


ISGRI 


PICsIT 


SPI 


39 


52677.2 - 52679.8 


149.3 


142.3 




- 


40 


52680.2 - 52681.6 


122.6 


56.0 


106.5 




41 


52683.2 - 52685.8 






204.0 




42 


52686.4 - 52688.2 


75.4 


61.6 




138.5 


43 


52689.6 - 52691.7 




30.7 




154.9 


44 


52692.2 - 52694.4 




9.8 




146.9 


45 


52695.2 - 52696.7 




2.1 




27.4 


102 


52866.3 - 52868.1 




37.9 




78.0 


103 


52868.6 - 52868.8 




21.6 






Total 




347.3 


362.0 


310.5 


545.7 



on board INTEGRAL. We investigate the presence of the 
440 keV line and thanks to the wide energy range cov- 
ered by INTEGRAL instruments and to the good statis- 
tics achieved by the summed data sets, we are able to per- 
form a detailed phase resolved spectroscopy on the Crab 
pulsed emission over the wide energy range 3-500 keV. 

2. Observation and Data reduction 

The Internatio nal Gamma-Ray Ast rophysics Laboratory 
(INTEGRAL; IWinkler et all bo03|) observed the Crab 
nebula and pulsar for calibration purposes several times 
from February 2003 (rev. 39) to August 2003 (rev. 103). 

The INTEGRAL payloa d is composed of three high 
energy instruments. JEM-X ijLund et al.ll2003|) consists of 
two identical coded-aperture mask telescopes with a geo- 
metrical area of 500 cm 2 and an angular resolution of 3 
arcmin across an effective field of view of about 10°. The 
detector at the focal plane, a Microstrip Gas Chamber, 
operates in the energy range 3-35 keV with an energy r es- 
olution of -17% at 6 keV. SPI l|Vedrenne et al.ll2003[) is 
a high spectral resolution gamma-ray telescope that con- 
sists of an array of 19 closely packed germanium detectors 
surrounded by an active anticoincidence shield of BGO. 
The imaging capabilities of the instrument are obtained 
with a tungsten coded aperture mask adopting a particu- 
lar observing strategy (dithering) . The fully coded field-of- 
view is 16°, and the angular resolution is 2.5°. The energy 
range extends from 20 keV to 8 MeV with a typical energy 
resolu tion of 2.5 keV at 1.3 MeV. IBIS (|Ubertini et all 
2003) is a coded apertu re telescope compos ed by two de- 
tection layers: ISGRI (|Lebrun et all l2003[) and PICsIT 
l|Di Cocco et alll2003|) . ISGRI is a large CdTe gamma-ray 
camera operating in the range 15 keV-1 MeV, with a geo- 
metrical area of 2621 cm 2 and an energy resolution of —8% 
at 60 keV. PICsIT is composed by 64 x 64 Caesium-Iodide 
(Csl) scintillation pixels working in the energy intervals 
175 keV - 10 MeV. 

The limited telemetry budged of INTEGRAL degrades 
the event timing information on board. Taking into ac- 



count all possible uncertainties affecting the time accuracy 
(the On Board Time (OBT) accuracy, the orbit predic- 
tion etc.L iKuiper et all l|2003|) estimated that the result- 
ing time resolution is 90/xs for IBIS, 130/is for SPI and 
150/iS for JEM-X, about 30-40% worse than the nominal 
time resolution of each single instrument; moreover the 
INTEGRAL absolu te timing accuracy, as estimated by 
iKtiiner et alJ (2003) from Crab data, is about 40/is. The 
IBIS/PICsIT detector cannot be routinely configured in 
photon-by-photon mode due to the high telemetry bud- 
get requested for this operational mode. For timing stud- 
ies then, observers can select the spectral-timing mode 
in which the whole detector counts are accumulated on 
board in up to eight energy bands (default 4) and with 
an integration time in the range 0.97-500 ms (default 3.9 
ms). 

The analysis performed in this paper uses on axis JEM- 
X observations and SPI data relative to pointings less than 
6°. IBIS observations have a maximum off- axis angle of 1°. 
They are relative to different configuration of the instru- 
ment: in particular, data are accumulated on board with 
different rise-time selections. We then verified that spec- 
tra relative to different science windows have negligible 
differences in the energy range considered for the spectral 
analysis. IBIS/PICsIT observation intervals with time res- 
olution of 1 ms. have been considered and to improve the 
statistics of the light curve, data from rev. 0041 with an 
off-axis of 9.6° were also included. 

Table 1 summarizes the log of the observations used in 
this analysis together with the relative time exposures. 

To obtain spectra and light curves of the sources 
present in the field, INTEGRAL official software (OS A 1 ) 
wh ose algorithms are described in lOoldwurm et al l 12003) 



Gros et alJll2003l) for TBTS. ISkinner fc Connelll <T2003h . 
Stronsl l|200aV ror SPI, and I Wester gaarH et alJl|20rB 



and 
and 

for JEM-X, repeats a shadowgram deconvolution process 
several times, by selecting events in the energy and time 



Available at 

http : //isdc .unige . ch/ index . cgi?Sof t+download 
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Fig. 2. Spectral index vs. phase measured by JEM-X in 
the energy range 3-20 keV (top panel), IBIS/ISGRI in the 
energy range 20-500 keV (middle panel) and SPI in the 
energy range 40-200 keV (bottom panel). 



intervals of interest. However, when the source positions 
have already been determined (because a priori known or 
predetermined by a shadowgram deconvolution), it is al- 
ternatively possible, to select only the detector pixels fully 
illuminated by the source. This method simplifies the ac- 
cumulation of phase resolved spectra. The amount of illu- 
mination from a given source, normalized to the maximum 
illumination value is called Photon Illumination Fraction 
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Fig. 3. JEM-X, IBIS/ISGRI and SPI spectrum of the 
phase interval 0.99-1.00 fitted with a single power law 



(PIF) and is generated by the standard software for JEM- 
X and IBIS/ISGRI. For these instruments, we run the 
standard pipeline (OSA vers. 4.2) up to the "DEAD" level 
that includes the conversion from detector energy channels 
(PHA) to energy channels corrected for instrumental ef- 
fects (PI), the selection for the Good Time Intervals (GTI) 
and the correction for the instrument dead time and we 
selected events with PIF= 1. However, IBIS/ISGRI con- 
version PHA-to-PI has been performed through our own 
calibration file generated by one of the authors (A.S.). 
This file is based on on-ground and in-flight calibration 
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data and represents an improvement respect to the stan- 
dard one (see Appendix A). 

No PIF selection is possible for SPI; we then extracted 
the list files relative to the whole field of view running the 
standard pipeline (OSA vers. 4.2) up to the "COR" level 
that produces corrected events selected for the GTI. 

The response matrices used in the analysis of JEM-X 
and SPI data are provided by the standard software. The 
ISGRI response matrix has been generated with our own 
software to take into account the new PHA-to-PI calibra- 
tion file, (see Appendix A). 

No response matrix is available for spectral timing data 
with IBIS/PICsIT; the presently available matrices are in 
fact suitable only for spectral imaging data. Data from this 
detector have not been included in the spectral analysis. 

JEM-X spectral analysis was performed in the energy 
range 3-20 keV, IBIS/ISGRI and SPI spectra were fitted 
in the range 20-500 keV and 40-200 keV, respectively. 
Errors quoted in the paper are relative to 1 a confidence 
level for one interesting parameter. 



Power Law 



H 1 h 



Curved Power Law 



Fig. 4. Frequency histogram of the reduced x 2 for 256 dof 
obtained fitting simultaneously JEM-X2 and IBIS/ISGRI 
spectra with a single power law (top panel) and with the 
curved model of Eq. (1) 



of P and P in GRO format for each observation 
(our data set spans several months) were derived 
from Jodrell Bank Crab Pulsar Monthly Ephemcris 
(http://www.jb.man.ac.uk/) using contemporary radio 
ephemeris. 
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Fig. 5. Best fit parameters a and b measured fitting si- 
multaneously JEM-X2 and IBIS/ISGRI with the curved 
model of Eq. (1) vs. phase 

Phase histograms of the Crab pulsar were evaluated 
for each instrument and each observation using the period 
folding techn i que an d adding the various offsets quoted in 
IWalter et all (|2003|) to correct the time relation derived 
by the INTEGRAL Science Data Center. The resulting 
phase histograms in six energy bands from 3 keV to 360 
keV are shown in Fig. 1 in absolute phase with a phase 
resolution ranging from 0.01 (0.33 ms) to 0.03 (1.1 ms) 
according to the available statistics and to the instrument 
time resolution. 

The well-known double peaked structure is prominent in 
all the profiles with a high statistical significance and the 
known evolution of the Crab pulse profile with energy can 
be observed: the relative intensity of the first pulse respect 
to the second increases with energies together with the 
level of the bridge emission. 



3. Timing Analysis 

Arrival times were converted to the Solar System 
Barycentre with the DE200 ephemeris. The values 



4. Spectral Analysis 

Crab pulsar phase histograms with 100 phase bins for 
JEM-X1, JEM-X2 and IBIS/ISGRI and 50 phase bins for 
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SPI were generated for each energy channel of each instru- 
ment and organized in an energy-phase matrix. A detailed 
phase resolved spectral analysis in the range (—0.1, +0.46) 
was performed by selecting spectra in phase intervals 0.01 
wide in the two main peaks and 0.02 in the interpeak 
for JEM-X1, JEM-X2 and IBIS/ISGRI; phase intervals 
of double size (0.02 at the two peaks and 0.04 in the in- 
terpeak) were considered for SPI data. The nebular emis- 
sion and the instrumental background for each phase re- 
solved spectra were subtracted evaluating them from the 
off-pulse level in the pase interval (+0.60,-1-0.80) and con- 
sidering that the contribution of the pulsed emission in 
this phase interval is negligible (jWeisskopf et alJ l2004). 

Energy channels are uniformly rebinned in agreement 
with the response matrices and in order to have a mini- 
mum bin content of 20 counts. 

Spectra of each instrument were first modeled with a 
single power law; low energy absorption has been included 
in the JEM-X f its, fixing the absorbing column to the val- 
ues derived by IWeisskonf et al.1 l|2004h . Spectra from the 
two JEM-X units were fitted simultaneously introducing 
as free parameter a factor to take care of the instrumen- 
tal systematics. The best fit values of this intercalibration 
factor are in the range 0.8-1.06 with an average of 0.98. 

JEM-X reduced \ 2 are generally acceptable: they span 
the range 0.8 1.2 (d.o.f 266) with only four values out of 
42 above this range (between 1.29 and 1.33). 

IBIS/ISGRI fits gave values of reduced x 2 between 0.8 
and 1.8 (122 dof) with 16 values greater than 1.26 over 
42 and only one above 1.5. These higher values of reduced 
X 2 can be considered acceptable because they are due to 
local residuals at low energies as expected from the level 
of accuracy of the systematics in OS A 4.2 software 2 . The 
SPI values of the reduced \ 2 li e m the expected range 
0.7-1.4 (41 d.o.f.) with only two values at 1.6. The best 
fit spectral indices are shown in Fig. 2 vs. phase together 
with the light curves for the three instruments. As example 
in Fig. 3 (top panels) the JEM-X, IBIS/ISGRI and SPI 
spectra relative to the the first peak are shown together 
with the residuals respect to the simple power-law model 
(bottom panels). 

The same phase dependence is clearly apparent in each 
plot in Fig. 2: the first peak has the softest spectrum, 
whereas the hardest emission is produced in the interpeak. 
The statistical significance of the softening of the spectral 
index has been evaluated fitting the photon index in the 
leading edge of the first peak with a constant and with a 
line. Applying the F-test to the derived x 2 , a significance 
of 99.6% in the JEM-X2 energy range and 99.7% in the 
IBIS/ISGRI 20-5 00 keV can be i nferre d co nfirming the re- 
sults o btained bv lPravdo et al 1 Jl997h and lMassaro et all 
l|2000h . 

The presence of a line at 440 keV in the ISGRI spec- 
trum relative to the phase interval (0.27-0.47) has also 
been investigated. We find a 3 a upper limit of 1.4xl0 -3 



ph cm 2 s 1 consistent with t he presence of the line de- 
tected bv lMassaro et all l|l99ll) . 

Comparing JEM-X, SPI and IBIS/ISGRI results, we 
note that spectral indices are clearly increasing with 
energy over all the phase intervals, in agreement with 
BeppoSAX results l|Massaro et alJ l2000l) . It is already 
known that the spectral energy distribution of the Crab 
pulsed emission is continously steepening f rom the optical 
freque ncies to 7-rays. In the X-ray range, IMassaro et alJ 
(2000) showed that a suitable model is the curved power 
law with a continously steepening described by the follow- 
ing formula: 



F{E) = K E 



-(a+bLog(E)) 



CD 



see the report osa_sci_val_isgri-L0.pdf available at 
http : / / isdc .unige . ch/Sof t/download/osa/osa_doc/prod 



where a corresponds to the photon index at 1 keV and b 
measures the curvature of the spectral distribution. 

We fitted at first JEM-X2 and ISGRI spectra si- 
multaneously with a single power law introducing a 
normalization factor in the model to take into account 
the intercalibration systematics between the two instru- 
ments. In this wide band fits we considered only JEM-X2 
units that shows a better calibration compared to the 
BeppoSAX-MECS. We find, in fact that the discrepancies 
in the photon indices measured by the two instruments, 
in the common energy range 2-10 keV, in the same phase 
intervals is lower than 5%. SPI has not been included 
because of the lower statistics. The resulted reduced 
X 2 have values generally unacceptable for the expected 
distribution with 256 degree of freedom. The distribution 
of the reduced x 2 1S shown in the top panel of Fig. 4: 
values > 1.3 are r elative to the s p ectra in the two main 
peaks. Following IMassaro et alJ l)2000fl approach, we 
fitted then JEM-X2 and ISGRI spectra simultaneously 
with the curved model of Eq. (1). All fits gave acceptable 
X 2 , as shown in the bottom panel of Fig. 4 and the best fit 
values of intercalibration factors are compatible with the 
constant 0.92+0.01 over all phase intervals. The best fit 
values of the two parameters vs phase are shown in Fig. 
5. The bending parameter b is statistically compatible 
with a single value over all phase intervals as found by 
Massaro et al. (2001) over 3 wider intervals. The fit with 
a constant gave a value of 0.14+0.02, where the error 
represents the spread around the average. 



5. Discussion 

INTEGRAL observations of the Crab pulsar provided a 
high-statistical data set for the study of the spectral and 
phase distribution of the pulsed emission over a wide X- 
ray energy interval (3-500 keV). The main result from the 
analysis of these data is that the pulsed spectral distribu- 
tion can be accurately represented by a curved function. 
The values of the bending parameter b in different phase 
intervals are consistent within the errors with a constant. 
This is a first independent confirmation of BeppoSAX re- 



i his is a nrst independent confirmation 01 tseppob AA re- 
sults presented in IMassaro et al . (2000). Other results of 
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the timing and spectral analysis can be summarized in the 
following points: 

— the spectral distribution changes with phase: the pho- 
ton index softens towards PI, hardens in the Ip region 
and increases again in the second peak with a char- 
acteristic reverse S shape over the 3-500 keV energy 
range; 

— the photon index at the PI leading edge shows a 
significative increase both in the JEM-X2 and in 
the IBIS/ISGRI e nergy ra n ges co nfirm ing the results 
obtain ed by IPravdo et al.l 1119971) and IMassaro et alJ 
l|2000l) yet questioned bv I Vivekanandl l|2002jh 

— the analysis of IBIS/ISGRI spectrum relative to the 
phase interval (0.27-0.47) does not rul e out the pres- 
ence o f the 440 keV line detected by IMassaro et all 
lll99lh . 

The value of the bending parameter (0.14±0.02) is sim- 
ilar to the values obtained for other Crab-like pulsars 
dde Plaa et all 120031: ICusumano et al.ll200lt iMineo et all 
|2004|) strongly suggesting a common characteristic of these 
sources. 

A log-parabolic spectrum, can be interpreted in term of 
the physics of the particle acceleration. It can be obtained 
when the acce l eration d ecreases with the particle energy 
IMassaro et alJ ()2004alb|h In the case of the pulsar envi- 
ronment, this could result from several crossings of the 
magnetosphere gaps with a time of permanence inside the 
acceleration region that decrease with the energy of the 
particles. 

The Spectral Energy Distribution (SED) of the log- 
parabolic law has a maximum at the energy E p given by 



(2) 



Considering the variation with phase of the parameter a, 
the values of E p ranges between ~10 keV and ^130 keV, 
in first peak and in the interpeak phase intervals, respec- 
tively. A possible explanation of the phase variation of the 
maximum energy is that we are observing photons emit- 
ted at different levels of the magnetosphere. T he three- 
dimensional outer gap model IjCheng et all2 000) seems to 
provide a viable theoretical description. In the framework 
of this model, curvature gamma-rays are converted into 
electron-positron pairs by the interaction with the mag- 
netic field and X-ray photons are then radiated by these 
secondary particles as synchrotron emission with a typical 
photon energy E syn : 



E 



syn{ r ) — 2 



E, 



he B(r) sin/3(r) 



mc 



(3) 



where E e is the electron energy B{r) is the dipole magnetic 
field, sin f3(r) oc (r) 1 / 2 the pitch angle and r the hight from 
the star surface within the magnetosphere. Assuming that 
E e is a fraction of the curvature energy estimated as: 



(4) 



where 7 e (V) oc (r) -1 / 8 is the local Lorentz factor and 
s(r) oc (r) 1 / 2 is the curvature radius, the ratio between 
the two SED maxima can be related to the ratio of the 
height of the emission regions n and r 2 : 



Esynjrx) _ Ecurjr-j) B(ri) sinfffri) 
E syn {r 2 ) E cur (r 2 ) B(r 2 ) sin/3 (r 2 )' 



(5) 



Introducing the dependence from r of each variable we 
find the following simplified relation: 



E S yn{ri) 

E syn {r 2 ) 



2.7 



(6) 



Results from INTEGRAL spectral phase resolved analysis 
on the two average energies imply a ratio of the heights of 
emittin g regions of ~ 0.3-0 .4 in agreement with the plot in 
Fig. 9 of Chen g et al 1 l|2000h where the level of the emission 
region vs phase computed for the Crab pulsar is shown. 

The present analysis of the Crab X-ray pulsed emis- 
sion confirms that a wide energy band analysis is very 
important for the study and understanding of the SED of 
radio pulsars. Future mission sensitive to the 7 rays should 
include this source as a primary scientific target. 

Acknowledgements. TM is grateful to Enrico Massaro for his 
helpful suggestions and discussion on the paper. The authors 
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that greatly improved the scientific content of the paper. 

Appendix A: ISGRI energy correction and 
response matrix 

The ISGRI detection layer, that consists of an f28xI28 
array of independent CdTe detector pixels, suffer of a 
rather severe " Charge Loss Effect" , common to this kind 
of detectors. To take into account th ese effects, the solu- 
tion adopted by the calibration team l|Lebrun et alJ l2003) 
is to perform an energy correction as a function of the 
pulse rise-time using multiplicative coefficients stored in a 
"Look-Up" Table (LUT2). The version of the LUT2 dis- 
tributed with the OSA software is not yet optimized 3 , and 
it introduces artificial features in the spectra, the most 
relevant in the 80 keV region, that are compensated with 
an ad-hoc modifications of the ISGRI effective area 4 (see 
left panel of Fig.^jL). However, the intensity of the arti- 
ficial features strongly depend on the spectral shape and 
the analysis of sources with spectral shapes different from 
that of the Crab might be affected. 

A new LUT2 5 based on on-ground and in-flight cal- 
ibration data has been generated by one of the authors 
A.Segreto. The better energy correction is confirmed by 



3 See A. Segreto talk at the Internal INTEGRAL workshop 
available at 

http : //www . rssd. esa. int/Integral/workshops/ Jan2005/ 

4 see the report osa_sci_valJsgri-L0.pdf available at 
http: //isdc .unige . ch/Soft /download/ osa/osa_doc/prod 

5 the file is available from the web page 

http : //www. if c . inaf . it/~f errigno/ integral/ISGRI_alternative_IC 
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Fig. A.l. ISGRI effective area vs energy for the standard software (left panel) and for the A.S. version (right panel). 
Features in the range 30-40 keV are due to CdTe absorption edges. 



the fact that it is no more necessary to introduce ad-hoc 
wiggles in the effective area, as shown in the right panel 
of Fig.^L. Moreover, the new effective area gives a value 
of the Crab spectral index in better agreement with the 
one quoted in literature and measured by the other in- 
strument on-board INTEGRAL 6 (see also Table EJ). In 
Fig. fXJi , the residuals of the power law fit of the Crab 
spectrum obtained processing the ISGRI data with the 
standard OSA 4.2 response matrix (left panel) and with 
the A.S. LUT2 and effective area (right panel) are shown. 

Tests on this matrix have been performed analysing 
sources which are detected with good statistics up to 100- 
200 keV. The spectral parameter derived with the two ma- 
trices are plotted in Table EJ. together with the values 
quoted in literature. The matrix we adopted gives best fit 
values of the spectral parameters in agreement with the 
ones quoted in literature and x 2 values generally lower 
than the one derived from the standard OSA 4.2 matrix. 
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Fig. A. 2. Left Panel: Residuals of the power law fit of the Crab spectrum obtained processing ISGRI data with 
the standard LUT2 and effective area. Right Panel: Residuals of the power law fit of the Crab spectrum obtained 
processing ISGRI data with the A.S. LUT2 and effective area 
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Table A.l. Comparison of the best fit parameters obtained using the standard matrix (ISDC) and the new A.S. 
version. 



Parameter OS A 4.2 


A.S. Version 


Other experiment 


Crab nebula+pulsar; Power Law 






a 2.19±0.03* 


2.09±0.03 


2.1 


K (ph cm" 2 s _1 KeV" 1 ) 10.2±0.2* 
X 2 (d.o.f.) 5.6 (95) 


9.6±0.2 


9.7 


1.3 (95) 


Ref (1) 


GX 1+4; Comptonization continuum model (COMPTT in XSPEC) 
To (keV) 1.27 (frozen) 1.27 (frozen) 1.27±i;t 


Te(keV) 13.5±0.5 


12.6±0.4 


10.3 + g 


T 2.7±0.1 


3.3±0.2 


3.7l| 


K (ph cm" 2 s" 1 KeV- 1 ) (2.1±0.1) x 1(T 2 
X 2 (d.o.f.) 1.5 (58) 


(1.9±0.1)xl0~ 2 


2.i± 2 jxi(n 2 


1.1 (58) 


Ref (2) 


Cygnus XI; Cut-off power law 






a 1.66±0.04* 


1.45±0.04 


1.45±0.01 


E c (keV) 158±12* 


122±7 


162±9 


K (ph cm" 2 s" 1 KeV -1 ) 1.7±0.3* 
X 2 (d.o.f.) 3.2 (48) 


1.2±0.3 


0.91±0.04 


1.9 (48) 


Ref (3) 


Sco XI; Bremstralung + Power law 






KT (keV) 4.2±0.1 


4.3±0.1 


4.51±0.08 


Flux2Q-wkev (erg cm -2 s" 1 ) (5.5±0.3) x 10~ 9 


(6.5±0.4)xl0~ 9 


(7.4±0.6)xl0^ 9 


a 2.9±0.2 


2.8±0.2 


2.4±0.3 


Flux 20 -2ookev [erg cm" 2 s' 1 ] (0.8±0.4) x HT 9 
X 2 (d.o.f.) 1.3 (52) 


(0.7±0.4)xl0~ 9 
1.4 (52) 


(1.04±0.08)xl0~ 9 
Ref (4) 



* Errors are compu t ed int rodu cing systematic error s in order to reduc e th e \ 2 below 2.0 
m lToor fc Sewardl l|l974l h (2) iGallowavl l|200(tl : f3l lDove et all l)l99Slh f4'l lD'Amico et alJ (|200lfr 



